High statistics study of the / (980) resonance in 77 — ► -k + -k production 
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We report on a high statistics measurement of the cross section of the process 
77 — > 7r + 7r~ in the it + tt~ invariant mass range 0.8 GeV/c 2 < W < 1.5 GeV/c 2 with 
85.9 fbr 1 of data collected at \/s = 10.58 GeV and 10.52 GeV with the Belle detec- 
tor. A clear signal for the /o(980) resonance is observed. From a fit to the mass spec- 
trum, the mass, n + n~ and two-photon decay widths of the resonance are found to be 
985.6 tl'l (stat) ±{\l (syst) MeV/c 2 , 34.2 t\l i (stat) ±1J (syst) MeV, and 205 ±|| (stat) t{" (syst) 
eV, respectively. 

PACS numbers: 13.66.Bc, 14.40.Gx 



The nature of low mass (below 1 GeV/c 2 ) scalar 
mesons has been a puzzle for decades with little progress 
made on its understanding [1]. Among the low mass 
scalar mesons, the existence of the fo (980) and ao (980) 
mesons is experimentally well established. One of the 
key ingredients in understanding their nature is measure- 
ment of the two-photon production cross sections and in 
particular the two-photon widths extracted from them. 
According to a relativistic quark model calculation, as- 
suming the fo (980) meson to be a non-strange qq state, 
its two-photon width should be in the range 1.3 keV to 
1.8 keV [2}. However, a much smaller width is expected 
for an exotic state (0.2 - 0.6 keV for a KK molecule 
state) 0, or for an ss state (0.3 - 0.5 keV) 0. 

A B factory is one of the best laboratories for a de- 
tailed investigation of low mass scalar mesons through 
two-photon production, where overwhelming statistics 
can be obtained. Two-photon production of mesons 
has advantages over meson production in hadronic pro- 
cesses; the production rate can be reliably calculated 
from QED with T 77 as the only unknown parameter. 
In addition, a meson can be produced alone without ad- 
ditional hadronic debris, and the quantum numbers of 
the final state are restricted to states of charge conju- 
gation C = +1 with J = 1 forbidden (Landau- Yang's 
theorem Q). 

In the past, using 209 pb" 1 of 77 — * 7r + 7r~ data, 
Mark II observed a shoulder in the 1 GeV/c 2 mass re- 
gion, which was tentatively identified as the fo (980) res- 
onance Q. The reaction 77 — * 7r°7r° was analyzed us- 
ing 97 pb _1 of data taken with the Crystal Ball detec- 
tor [3] ■ They found a hint of fo (980) formation with 
a significance of 2.2 standard deviations. Measurements 



of 77 — > 7r°7r° were also performed with the JADE de- 
tector using 149 pb -1 data [{|. They observed a small 
shoulder at around 1 GeV/c 2 , which was interpreted as 
the production of the fo (980). CELLO studied the re- 
action 77 — ► 7r + 7r~ using a data sample of 86 pb _1 and 
concluded that an / (980) signal at the level reported in 
Refs. @, 0, [f| cannot be excluded with their errors 0] . 

Using data from Mark II, Crystal Ball, and CELLO, 
Boglione and Pennington (BP) performed an amplitude 
analysis of 77 — * ir + n~ and 77 — * ir°ir cross sec- 
tions [10(. They found two distinct classes of solutions 
where one solution has a peak ( "peak solution" ) and the 
other has a wiggle ("dip solution") in the fo (980) mass 
region. The two solutions give quite different results for 
the two-photon width of the fo (980) and the size of the 
S-w&ve component. Thus, it is important to distinguish 
them experimentally 

In this paper, we report on a high statistics study of the 
/o(980) meson in the 77 — > n + ir~ reaction based on data 
taken with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider [ill ]. The data sample corresponds 
to a total integrated luminosity of 85.9 fb _1 , accumulated 
on the T(45) resonance {*Js = 10.58 GeV) and 60 MeV 
below the resonance (8.6 fb^ 1 of the total). Since the 
difference in the cross sections between the two energies 
is only about 0.3%, we combine both samples. We ob- 
serve the two-photon process e + e _ — > e + e - 7r + 7r~ in the 
"zero-tag" mode, where neither the final-state electron 
nor positron is detected, and the 7r + 7r~ system has small 
transverse momentum. 

A comprehensive description of the Belle detector is 
given in Ref. [13]. Charged track coordinates near the 
collision point are measured by a silicon vertex detector 
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(SVD) that surrounds a 2 cm radius beryllium beam pipe. 
Track trajectory coordinates are reconstructed in a cen- 
tral drift chamber (CDC), and momentum measurements 
are made together with the SVD. An array of 1188 silica- 
aerogel Cherenkov counters (ACC) provides separation 
between kaons and pions for momenta above 1.2 GeV/c. 
The time-of-flight counter (TOF) system consists of a 
barrel-like arrangement of 128 plastic scintillation coun- 
ters and is effective for K/n separation for tracks with 
momenta below 1.2 GeV/c. Low energy kaons and pro- 
tons are also identified by specific ionization (dE/dx) 
measurements in the CDC. Photon detection and energy 
measurements of photons and electrons are provided by 
an electromagnetic calorimeter (ECL) consisting of an 
array of 8736 CsI(Tl) crystals all pointing toward the in- 
teraction point. These detector components are located 
in a uniform magnetic field of 1.5 T provided by a su- 
perconducting solenoid coil. An iron flux-return located 
outside the solenoid coil is instrumented to detect K\ 
mesons and to identify muons (KLM). 

Signal candidates are primarily triggered by a two- 
track trigger that requires two CDC tracks with asso- 
ciated TOF hits and ECL clusters with an opening angle 
greater than 135 degrees. Exclusive e + e~ — > e + e~7r + 7r~ 
events are selected by requiring two oppositely charged 
tracks coming from the interaction region; each track is 
required to satisfy dr < 0.1 cm and \dz\ < 2 cm, where 
dr (dz) is the r (z) component of the closest approach 
to the nominal collision point. The z axis of the de- 
tector is defined to be opposite to the direction of the 
positron beam and r is the transverse distance from the 
z axis. The difference of the dz's of the two tracks must 
satisfy the requirement \dz+ — dz—\ < 1 cm. The event 
must contain one and only one positively charged track 
that satisfies p t > 0.3 GeV/c and -0.47 < cos6» < 0.82, 
where pt and 9 are the transverse component of momen- 
tum and the angle with respect to the z-axis. The scalar 
sum of track momenta in each event is required to be 
smaller than 6 GeV/c, and the sum of the ECL energies 
of the event must be less than 6 GeV. Events should 
not include an extra track with p t > 0.1 GeV/c. The 
cosine of the opening angle of the tracks must be greater 
than —0.997 to reject cosmic-ray events. The sum of the 
transverse momentum vectors of the two tracks (X^Pt) 
should satisfy | J^Pt l< 0-1 GeV/c; this requirement 
separates exclusive two-track events from quasi-real two- 
photon collisions. 

Electrons and positrons are distinguished from hadrons 
using the ratio E/p, where E is the energy measured 
in the ECL, and p is the momentum from the CDC. 
Kaon (proton) candidates are identified using normal- 
ized kaon (proton) and pion likelihood functions obtained 
from the particle identification system {Lk (L p ) and L„, 
respectively) with the criterion Lk/{Lk + L n ) > 0.25 
(Lp/(L p + L^) > 0.5), which gives a typical identifica- 
tion efficiency of 90% with a pion misidentification prob- 
ability of 3%. All charged tracks that are not identified 
as electrons, kaons or protons are treated as pions. We 



require both tracks to be pions. 

In this measurement, the KLM detector cannot be 
used for muon identification, since it is insensitive in 
the region of interest where the transverse momenta 
of tracks are below 0.8 GeV/c. Therefore, we have 
developed a method for statistically separating ir + ir~ 
and n + pT events using ECL information; muons de- 
posit energy corresponding to the ionization loss for 
minimum ionizing particles, while pions give wider en- 
ergy distributions since they interact hadronically in the 
ECL, which corresponds to approximately one interac- 
tion length of material. Probability density functions 
(PDFs) for the distributions of energy deposits from 

M + M " (tt+tt-) pairs P« M _ (E + ,E_) (P« ff _ (£+,£_)) 
are obtained with GEANT-3 [jj] Monte Carlo (MC) sim- 
ulation. Here i represents the i-th bin of (W, |cos#*|) in 
20 MeV/c 2 and 0.1 steps, where W is the invariant mass 
of the 7r + 7r~ (or /i + /i~) pair in each event (the pion mass 
is assumed in the calculation), and 9* is the polar angle 
of the produced tt^ meson (or /x^ lepton) in the center- 
of-mass system of two initial photons. Note that the 
effect of muons from pion decays is taken into account 
by the pion PDFs using this method. We obtain A l \ the 
fraction of in the i-th bin through the equation: 

+(l-r«)P^V . 

where N^ ata {E + , E-) is the distribution of data and N^ t 
is the total number of events in that bin. The values of ra- 
tios rW obtained must be corrected since the MC cannot 
simulate hadronic interactions accurately enough. By in- 
troducing mis-ID probabilities, P-m,-^^ and P^-t^, the 
r value for each bin (the bin number i is omitted) can be 
written as: 

N ~t~ A^ 7r7r P 7r7r — >^ -^Lt/^P/^jU — > 7T7T 

AW + N ™ ' 

where (N^) is the number of true (7r + 7r~) 
pair events in that bin. We assume that P?™— >^4 and 
P////-nrff are independent of W. Applying the h/tv sepa- 
ration method described above to a sample of data events 
positively identified as muons by the KLM, we find that 
P^-mir is statistically consistent with zero. The values 
of P^-^^ in each | cos 9* | bin are determined such that 
the ratio of the data and MC for fi + fi~ pairs, which is one 
ideally, gives a straight line in the W spectrum. The val- 
ues of Ptttt-*^ vary between 0.08 to 0.13 in | cos6**| bins. 
Because they are determined for each bin of | cos#*|, the 
bin-by-bin variation of systematic errors is rather large 
in the angular distribution. 

The total cross section for 77 — > ir + ir~ with | cos6**| < 
0.6 is evaluated using the following equation: 

a ^ + *- = etl . g -e dct -AW-^-jLdt- 
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Here AN e + e -^ >e + e -„+ 7r - is the number of events in a W 
bin, -j^r is the two-photon luminosity function [TjJ and 



/ Ldt = 85.9 fb~ is the integrated luminosity. The size 
of the W bin is chosen to be 5 MeV/c 2 ; a typical mass 
resolution for a ir + ir~ system is 2 MeV/c 2 according to a 
MC study. The detection (trigger) efficiencies, edet (etrg) 
are estimated with a MC simulation. Events of the pro- 
cess 77 — > 7r + 7r - are generated using TREPS [15]. The 
detection efficiency is extracted from MC simulation and 
the trigger efficiency is estimated with the trigger simula- 
tor. Since the trigger simulator does not simulate triggers 
well, particularly in the low energy region, the efficiency 
values have to be corrected. We calculate the correction 
factors by comparing e + e~ — > e + e~e + e~ events in data 
and MC that are triggered by the two-track trigger. The 
resulting factors steeply rise from 0.5 at W — 0.8 GeV/c 2 
to 0.8 at W = 1 GeV/c 2 and then increase gradually for 
higher W . The muon-background subtraction and all the 
correction factors are applied using smooth functions ob- 
tained by parameterizing the results of bin-by-bin anal- 
yses. 

The total cross section obtained is shown in Fig. [1] to- 
gether with the results of some past experiments; an ex- 
panded view of the /o(98Q) region is shown in Fig. [2][a). 
A clear peak corresponding to the /o(980) meson is vis- 
ible, and thus the peak solution of the BP analysis is 
selected. Systematic errors for the total cross section are 
summarized in Table U They are dominated by the un- 
certainty of the fi/ir separation and that of the trigger 
efficiency. Systematic errors arising from the fi/ir sep- 
aration are estimated by changing the value P^-^, in 
the allowable range in each angular bin. Since 
events are well identified for W > 1.6 GeV/c 2 , the al- 
lowable range is determined in this region. These well 
identified events are also used in estimating sys- 

tematic errors of the trigger efficiency. Comparing data 
and MC for fi + n~ events in the region W > 1.6 GeV/c 2 
and extrapolating linearly downward, the systematic er- 
rors are found to be 4% at W = 1.5 GeV/c 2 and 10% at 
W = 0.8 GeV/c 2 . The total systematic error is obtained 
by summing the systematic errors in quadrature and is 
also shown in Fig. [TJ 



TABLE I: Summary of systematic errors for the 77 — > ir + ir~ 
cross section. A range is shown when the uncertainty has W 
dependence. 



Parameter 


Syst. error (%) 


Tracking efficiency 


2.4 


Trigger efficiency 


4-10 


K 1 7r-separation 


- 1 


/i/7r-separation 


5-7 


Luminosity function 


5 


Integrated luminosity 


1.4 


Total 


11.1 - 12.3 



+ Belle 
--■ syst. errors 
□ Mark II 
A CELLO 



I cost?" I < 0.6 



8 9 1 1.1 12 1 3 1 4 1,5 

W (GeV/c 2 ) 

FIG. 1: The total cross section of e + e~ — > e + e - 7r + 7r~ be- 
tween 0.8 and 1.5 GeV/c 2 for cos 6>* < 0.6. The Belle 
data are represented by crosses with statistical error bars, 
the CELLO data are the triangles and the Mark II data are 
squares. Dashed lines are upper and lower systematic uncer- 
tainties for the Belle data. 



Our results are in good agreement with past experiments 
except for the fa (1270) mass peak region, where they are 
about 10 to 15% larger, but still within the systematic 
errors. 

A fit to the 77 — > 7r + 7r~ total cross section is performed 
to obtain the parameters of the /o (980) meson. We have 
to take into account the effect of the KK channel that 
opens within the fo (980) mass region. The fitting func- 
tion for the scalar resonance /o(980) is parameterized as 
follows: 



F fo t 



r BG 



r BG 



(1) 



where the factor 4.8 includes the fiducial angular accep- 
tance |cos6>*| < 0.6, (3 X = y 1 - ^wr- is the velocity 
of the particles X with mass Mx in the two-body fi- 
nal states, is the amplitude of the /o(980) meson, 
which interferes with the helicity-0-background ampli- 
tude VOf) w hh relative phase (/?, and cr BG is the to- 
tal background cross section. The amplitude T'° can be 
written as 

77-' (2) 



16tt 



where gj xx is related to the partial width of the /o(980) 
meson via Txx(fo) = ~Tg§Mj^ ■ The factor Dj Q is given 
as follows 16]: 



D f0 (W) = M% 



w 2 + (M fo )-nl° (w) 
5ffl£ (M /o )-n£ (w) , 
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where for X = n or K . 



n£°(w) = 



Px9f xx 



16tt 



I A In 

7T 1+Px 



(3) 



The factor /3k is real in the region W > 2Af^ and 
becomes imaginary for VF < _2Mk- The mass differ- 
ence between fC ± and if (if ) is included by using 
Pk = \(Pk± + Pk )- 

In the fit, we assume <Tq to be constant and the rel- 
ative phase to be a slowly varying function of W; this is 
motivated by the nearly energy- independent behavior of 
the scalar Born amplitude [17|. We fix g 2 gKK /5/ o;T7r = 
4.21±0.25 (stat)±0.21 (syst) taking the latest value from 
the BES measurement [18J. The background function 
(T BG is evaluated by fitting the cross section with a 4- 
th order polynomial in W outside of the fo (980) region 
0.85 GeV/c 2 < W < 0.93 GeV/c 2 and 1.03 GeV/c 2 < 
IF < 1.15 GeV/c 2 . The value of x 2 /ndf for the fit is 
0.88 for 46 degrees of freedom (ndf). A fit to the fo (980) 
resonance is then performed with Eq. (fTJ) in the region 
0.93 GeV/c 2 < W < 1.03 GeV/c 2 , where the parameters 
of cr are fixed; the free parameters are Mf , <7/ o7r 7r, T 77 
(evaluated at the / (980) mass), a^ G and ip. 

The result of the fit is shown in Figs. Ela) and[2jb). 
In Fig. HJb), one can see a significant interference effect, 
which is visible as a deviation from a Breit-Wigner-like 
shape in Fig. [2fa). In the same figure, the cross sec- 
tion er(77 — > /o(980) — > Jf + Jf~) is also plotted, which 
is obtained by evaluating the first term in Eq. (fT|), sub- 
stituting (3k instead of and in Eq. ^ 9f KK instead 



of 



ff/o7rn 



Note that the cross section is zero below the 



threshold even though the amplitude is non zero. The 
value of x 2 /ndf of the fit is 1.04 for 15 ndf . The helicity- 
0-background component that interferes with the /o(980) 
meson (<Tq G ) is found to be 3.7I^; 2 nb. The value of ip is 
approximately 7r/2, which is consistent with the general 
phase shift study [10| . 

The parameters of the fo (980) meson are found to be 



r+rr- (fo) 
r 77 (fo) 



985.6 
34.2 
205 



t\l (stat) t\i (syst) McV/c 2 



±i?:l (stat) ±|;| (syst) MeV 



+95 (stat) i^(syst) eV. 
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The two-photon width given by the PDG [19( is 
T 77 (fo) — 310 (stat) eV, and the value found by 
BP is 280t?° eV. Our results are consistent with them 
within errors as well as with the prediction of the four- 
quark model of 270 eV 

The dominant systematic errors come from fitting. 
The value of T 77 (fo) is quite sensitive to changes in pa- 
rameters of the background cross section (fitted outside of 
the /o(980) resonance). Systematic errors are evaluated 
by changing each background parameter by ±ler, taking 
their correlations into account; the error is strongly cor- 
related with that of gf o7r -K (i.e. r^-t^- (fo)). The error 
in the normalization of the total cross section has lit- 
tle effect on the value of the /o(980) mass, however it 
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FIG. 2: Fitted curve: (a) shows the total cross section and 
(b) shows contributions of the resonance ((7(77 — > /o(980) — > 
7T + 7r - )) (solid line) and the interference (dashed). The cross 
section of (7(77 — > /o(980) — ► K + K~) is also shown (dotted). 



is a significant contribution to the error in T 77 (fo) and 
(fo)- The errors in gj oKK /gj o7r7r are also taken 
into account in the systematic errors. Individual system- 
atic errors are summed in quadrature to obtain the total 
uncertainty. 

In summary, we have made a high statistics measure- 
ment of the 77 — > 7r + 7r~ cross section in the 7r + 7r~ in- 
variant mass region 0.80 GeV/c 2 < W < 1.5 GeV/c 2 
in fine bins of W (5 MeV) and cos6»* (0.05) with the 
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Belle detector at the KEKB accelerator. We have ob- 
served a significant signal corresponding to the f (980) 
resonance. Our data clearly select the peak solution of 
the Boglione-Pennington amplitude analysis The 
total cross section is fitted to obtain the parameters of 
the /o(980) meson. Its two-photon width is found to be 
205 (stat) 1^7 (syst) eV, consistent with past exper- 
iments. 
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